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ABSTRACT 



Deep Very Large Array imaging of the binary X-ray source SS 433, sometimes 
classified as a microquasar, has been used to study the intrinsic brightness dis- 
tribution and evolution of its radio jets. The intrinsic brightness of the jets as 
a function of age at emission of the jet material r is recovered by removal of 
the Doppler boosting and projection effects. We find that intrinsically the two 
jets are remarkably similar when compared for equal r, and that they are best 
described by Doppler boosting of the form as expected for continuous jets. 

The intrinsic brightnesses of the jets as functions of age behave in complex ways. 
In the age range 60 < r < 150 days, the jet decays are best represented by 
exponential functions of r, but linear or power law functions are not statistically 
excluded. This is followed by a region out to r ~ 250 days during which the 
intrinsic brightness is essentially constant. At later times the jet decay can be fit 
roughly as exponential or power law functions of r. 

Subject headings: binaries: close — radio continuum: stars — stars: individual 
(SS433) 
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INTRODUCTION 



The galactic X-ray binary source SS433, consisting of a stellar-rnass b 



2008 



Bowler 



orbit about an ear ly-type star (IBlundell. Bowler. &: Schmidtobreick 



2008 



ack hole in close 



Hii: 



wig fc Gies 



1999), and is 



20101), is a miniature analogue of an AGN ( iMirabel &: Rodriguez 
often classified as a microquasar. Two mildly relativistic jets emerge from opposite sides of 
the compact object at speeds of 0.26 c. Modeling of the optical spectrum s hows that the jet 



system precesses wit 


1 a period of 162 days about a cone 


1979; 


Fabian & Rees 


1979; 


Milerom 


1979 


Mareon 


1984) 



1984J ). Imaging by the Very Large Array 



(VLA) conf irmed this picture, as the radio images sho wed helical jets on both sides of 



the source fIS pencer 



1979 



Hjellming fc Johnston 



1981a 



b|). Higher resolution 



by VLBI reveal the structure down to a scale of a few AU (jVermeulen et al. 



imag e s mad e 



19871 . 



19931 ). 



Analysis of the 15 GHz VLA-s cale structure of the j ets in SS433, with angular resolution 



of about 0.1", was presented in 



Roberts et al. 



(|2008[ ). Mult i-epoch dual-frequency analy sis 



of SS 433 during the summer of 2003 will be presented in 



( 120101 ). hereafter Paper HI, and for the summer of 2007 in 



R,oberts. Wardle. fc Ma! 



Bell. Wardle. fc Roberts 



lory 



mm, 



hereafter Paper IV. 

In this paper we use high dynamic range VLA images of SS 433 to study the radiative 
intensity of the two jets as a function of the material's birth epoch t and age at emission r 
(hereafter simply "age"); see Appendix 1 ([|7j) for our definitions of these quantities. Our 
goals are (i) to determine if the two jets are intrinsically the same, and (ii) to learn if the 
jets behave as individual non-interacting components or as a continuous stream. SS433 
offers a unique opportunity to answer these questions because it presents two jets with 
ever-changing mildly-relativistic velocities known as functions of time and position on the 
sky from their optical properties. In ^ we describe the observations and data reduction. In 
^we determine the properties of the jets and we discuss the physical implications of these 
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results in §H Our conclusions are summarized in ^ 



OBSERVATIONS AND DATA REDUCTION 



Interferometer data for 2003 July 11 (JD 2452832) were obtained from the VLA data 
archive. The array was in the A configuration with 27 working antennas. The frequency 
used was 4.86 GHz, with a bandwidth of 50 MH z for each of the two IF systems. The 



3 



2004J ) using the nearby source 



data were edited and phase calibrated in AIPS (INRAO 
1922+155, and amplitude calibrated against 3C48; the data were then imaged using the 
tasks IMAGR and CALIB, utilizing many cycles of pha s e and amplitude self-calibration. 



These data were p reviously used by 



of the jets, and by 



Miller- Jones et al. 



undell fc Bowlerl (12004J ) to study the symmetries 



( 120081 ) to study the magnetic field configuration in 



the jets. Figs. [U^a) & (b) show the distribution of total intensity for SS433, made with 
uniform weighting (ROBUST = —5); this image has resolution of 0.32". The kinematic 
model without nutation or velocity variations is shown in (a), and the model with nutation 
and velocity variations is shown in (b). The kinematic model used to define the jet locus 
is described in Appendix 1, ^ A naturally weighted image with a resolution of 0.47" 
(ROBUST = +5), shown in Figure [21 reveals very low level emission out to distances of at 
least 6", corresponding to ages of about 800 days for both jets. 

For all of the analysis that follows we used the image shown FiglH and included the 
velocity variations and nutation. This accounts for the "spiking" apparent in many of the 
figures. Uncertainty in the position of the model jet locus is by far the largest source of 
fiuctuations in the total intensity curves; the root-mean-square difference between profiles 
generated with and without these terms is about 1.5 mJy/beam. Thermal noise is negligible 
until ly falls below about 0.1 mJy/beam. 
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Comparison of the images with the kinematic model shows them to be compatible, with 
the kinematic locus being the leading edge of the jets; the source also contains significant 
off-jet material. We see no unambiguous evidence of either nutation or significant jet 
velocity variations, but this is not surprising given the limited resolution of the images. 

3. PROPERTIES OF THE SS 433 JETS 

3.1. Normalization Factors 

In order to study the intrinsic properties of the jets, we used the kinematic model 
to determine the effects on the observed jets of projection and Doppler beaming. In 
what follows, we will present jet properties as functions of age r instead of the birth 
epoch t because we expect the aging of the jet material to be the most important factor 
determining its properties. Figure |3] shows r for both east and west jets as functions of t, 
and it can be used to estimate r at any location on the images, using the average proper 
motion of ~ 8 mas d~^. Projection effects were determined using model jets consisting 
of discrete components equally-spaced in birth epoch t; the projected density on the sky 
was determined by performing a beam-averaged count of the number of components per 
beam area as functions of position down the kinematic locus of the the two jets. This 
was normalized to unity at the core. The Doppler factor D of each component was found 
from the kinematic model, and Doppler boosts calculated as D°'~^"', where a continuous 
jet has n = 2 and an isolated component has n = 3 (see Appendix 2 ([|8])), and a ~ 0.7 
is the spectral index of the jet {F^, oc Paper III, Paper IV). The total normalization 
factors were obtained by performing a beam-weighted sum of the boosts as a function of 
position down the model jets; all quantities were determined as functions of r. Figure H] 
shows combined normalizations for the jets in SS 433, expressed as multiplicative factors to 
be applied to observed total intensities, as functions of r, for n = 2. Curves for n = 3 are 
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similar and not shown. Data from Figure H] will be used below to create model jet intensity 
curves. The ratio of the normalizations for the two values of n ranges between 0.92 and 
1.12, sufficiently different from unity to distinguish the two in the regions r < 150 days 
where the jet flux is greatest. 

3.2. Determination of the Intrinsic Structure 

Figure E] shows the profiles down the model locus of the observed total intensities of 
the jets derived from the image in Figure [1] and the naturally-weighted image in Figure 121 
as functions of r. For the region r > 60 days, contamination by the core is less than 1% of 
the jet intensity in Figure [H the beam occupies the first 75 days in Figure El The loops in 
the curves for both jets are the result of the fact that material of multiple different birth 
epochs can have the same age as a result of differing light travel times. 

It is useful to know that the places along each jet where the radial velocity 
switches sign and the jet motions are both in the plane of the sky are located at 
t = T = 14, 125, 176, 288, 338, 450, and 500 days. At these ages the projection factors and 
Doppler boots of the two jets are equal, so if the two jets are intrinsically the same, the 
raw total intensities should be the same. Examination of the intensity curves shows that, 
allowing for the beams of full-width at half-maximum about Ar ^ 40 days for Figure [5]^a) 
and ~ 60 days for Figure El^b), this is satisfied. 

As can be seen from Figure [6l neither normalization for projection effects alone nor 
for Doppler boosting alone results in east and west jets appearing the same. This is true 
for either choice of n. Figure [7] shows the intrinsic brightnesJl] profiles of the jets in SS 433 

^We use this term to refer to the total intensities that we would observe were we in a 
frame co-moving with a piece of the jet. 
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created by normalizing for both projection effects and Doppler boosting for n = 2 and 
n = 3, as functions of r, derived from the image in Figure [H Especially in the case n = 2, 
the intrinsic brightnesses of the two jets are very similar when compared for equal ages. 
The small discrepancies between the derived intrinsic brightnesses of the east and west 
jets for n = 2 in the range 80 < r < 220 days can be understood as artifacts due to the 
tight bend in the east jet at RA ~ 1 and the first loop in the west jet, RA ^ — 1.2 ; these 
correspond to ages 80 < r < 130 days and 140 < r < 220 days, respectively. In the first 
range the measured total intensity for the east jet is artificially high, as seen in Fig. El^a), 
due to two parts of the jet locus being "in the beam" at the same time. In the second 
range, the total intensity of the west jet is falsely measured to be high, again corresponding 
to Fig. [7](a). The normalization process cannot account for the contamination of one side 
of the loop by the part of the other that is in the same beam (the procedure lacks a priori 
information about the ratio of the two contributions). This is minimal at the cusp of the 
loop, r ~ 180 days. However, when we reconstruct the west jet from the intrinsic brightness 
of the east jet and the normalization model this is accounted for, and the artifacts largely 
disappear (see below). While the same arguments apply to n = 3, they cannot account 
for the discrepancy over 60 < r < 90 days; the agreement for n = 3 is simply not as good 
where the jet is brightest, suggesting that n = 2. 

We can see the results of the normalizations in another way if we examine the total 
intensity ratio of the two jets. Figure |8] shows the ratio of observed total intensities in the 
form east jet divided by west jet and the predicted ratios for complete correction for n = 2 
and n = 3, all as functions of r. The fit of the combined normalization to the data is quite 
good for n = 2; the prediction for n = 3 is not as good a fit. In fact, the discrepancies 
between the data ratio and the n = 2 model ratio are explained in the same way described 
above. In the range 80 < r < 130 days, the east jet is measured to be artificially high, 
raising the measured ratio above the model ratio, in agreement with the figure; in the range 
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140 < r < 220 days, the west jet is measured artificially high, producing the opposite effect, 
again in accord with the figure. 

We can test the similarity of the two jets in a third way, by using the derived intrinsic 
brightness of the east jet in Figure [7| and the model normalization factors in Figure H] to 
predict the observed total intensity of the west jet. Comparisons for n = 2 and n = 3 are 
shown in Figure [9] for ages out to 350 days. In the region r < 350 days the east jet is a 
single-valued function of age, so it can be modeled (normalized) uniquely, and the loops in 
the west jet reconstructed free of artifacts. The small uncertainties in the intrinsic brightness 
of the east jet described above produce small disagreements over 60 < r < 200 days, as 
expected. The overall fit is better for n = 2 than for n = 3, also in agreement with the 
results above. Because t becomes a multiply-valued function of r beyond ages of about 
r = 350 days, it is not possible to reconstruct the west jet uniquely beyond this point. 

Finally, if we regard n simply as a phenomenological parameter, examination of 
Figs. U\M shows that n = 2 is favored over n = 3 because agreement of the intrinsic 
brightness of the two jets where they are brightest is superior. 

In summary, we conclude that in SS 433 jets the observed total intensities are compatible 
with the twin jets being intrinsically identical, and behaving as continuous jets. We address 
the question of the expected behavior of the observed total intensity on the Doppler boosts 
in Appendix 2 (§8]). 

4. DISCUSSION 

Three interesting questions about the jets are: (i) do they contain features, such as 
breaks or peaks in the intensity profiles in either or both jets that are not the effect of 
projection and Doppler boosting on an otherwise smooth intensity distribution, (ii) if so. 
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are they the same in the two jets, and (iii) if so, how do they behave in time. Such features 
might arise at the ejection of the jets (the core after all is quite variable; Paper III), or as 
functions of their aging and/or propagation. Figure [HJ] shows the intrinsic brightness for 
n = 2 and function of birth epoch. Because this compares east and west jet 

material of different ages there is no a priori reason to expect matching features. However, 
the regions of roughly constant brightness located over the range 150 < t < 250 days in 
both jets may be such features because we see the same behavior over the entire summers 
of 2003 (Paper III) and 2007 (Paper IV), but at different age ranges. Thus the flattening 
of the decay curve may be the result of variations in the central engine. Comparison of the 
two jets is hindered by our inability to extract unique properties of the west jet through the 
first tight loop. Nonetheless, what does seem clear is that after the initial rapid falloff of 
each jet, there is a leveling-off, and then a somewhat slower decline with large observational 
uncertainties. 

Assuming that the jets are dominated by aging rather than by variability in ejection 
properties, we can parameterize the physical processes that govern these different time 
ranges by fitting the normalized intensity curves with linear, exponential, and power law 
models. Figure [11] shows fits to the intrinsic brightness B of the jets for n = 2 over the 
range 60 < r < 150 days. Exponential fits of the form B oc e""^/^ yield half-lives of 
Ti/2 = Tin 2 = 41 days and 39 days for the east and west jets, respectively. Linear fits, 
while unrealistic for an entire jet, give half-lives of T = —B/2B of 46 days and 40 days, and 
power law fits yield exponents of —1.7 and —1.8, for the east and west jets, respectively. 
Reassuringly, these timescales are very similar to the results of multi-epoch imaging that 
follows the aging of specific pieces of the jet (Paper III, Paper IV). The root-mean-square 
deviations of the exponential fits of B{t) for the east jet and west jet are 1.0 mJy/beam and 
0.8 mJy/beam, which suggests that a reasonable lower limit to the uncertainties in the total 
intensity curves is about 1 mJy/beam, comparable to the difference between using model 
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curves with and without nutation and velocity variations (Q. Given the hmited span of 
these data, we are unable to distinguish definitively among these models; the exponentials 
appear superior, but this is not statistically significant. 

For constant speed, as in SS 433, models that predict intensity behavior with distance 
should be compared to data fit in age. The model of a freely-expanding spherical cloud 



of magnetized relativistic plasma (Ivan der Laan 



19661 ). which might be appropriate if 



the pieces of the jet behave as discrete non-interacting components, predicts that in the 
optically-thin regime the resulting synchrotron total intensity falls off as r~^^, where p is 
the electron energy distribution exponent. Since here p ~ 2.4, this predicts a decay with a 
power law in component radius of index ~ —4.8. Were the components freely expanding we 
w ould have r oc r, which would then be ruled out by the data. In the conical jet model 



of iHjellming fc Johnston] (1l988l ). the intensity of the jet falls off as dl oc z~"^dz, where 
z is distance down the jet, m = {7p — l)/(6 + 66), and 6 = and 6=1 correspond to 
freely-expanding and slowed expansion cases, respectively. For p = 2.4, 6 = corresponds 
to m = 2.6, and 6 = 1 to m = 1.3. Power fits to our data in this age range lie between the 
two cases. We also fit exponentials and power laws to the data for 250 < r < 800 days, and 
find for both jets T1/2 — 80 days or a < 4 where B oc r^'^, marginally incompatible with the 
freely expanding sphere model. 

5. CONCLUSIONS 

The principal results in this paper are: 

1. We have used a deep VLA A-array image of SS433 at 4.86 GHz to study the intrinsic 
brightness profiles of the twin jets. 

2. Radiation from both jets is detected out to at least 6" from the core, corresponding 
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to jet ages of about 800 days. 

3. The observed brightnesses of the jets are strongly affected by projection effects and 
Doppler boosting. 

4. Intrinsically the two jets are remarkably similar, and they are best described by 
Doppler boosting of the form as expected for a continuous jet. 

5. The intrinsic brightness of the jets behaves in a complex way that is not well described 
by single linear, exponential, or power law decay. 

6. During their first ~150 days, the jet decays are well represented by linear or 
exponential functions of age, with linear half-lives or exponential half-lives of about 
40 days, the same for the two jets. Power law fits to the data in this age range give 
exponents of about —1.8. 

7. There is a transition region, corresponding to jet ages between about 150 and 

250 days, during which the jets maintain roughly constant intrinsic brightnesses. This 

represents nearly one complete precession period. This also corresponds to about 
150 <t< 250 days in either jet. 

8. At later times the jet decay can be roughly fit as exponential functions of age, with 
exponential half-lives of about 80 days, or as power laws with indices of a < 4. 
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7. APPENDIX 1 



We used the geometric model o: 



parameters of 



Eikenberry et al. 



(1200 ll ). and the distance of 



( 120071 ). augmented with iet ye . 



nutation following 



Katz et al 



ocity variations following 



' iHjellming fc JohnstonI (|l981br). the precession 



Lockman. Blundell. fc Goss 



Blundell fc Bowlei] ( 120051 ) and 



( 119821 ). to determine the locus and kinematic properties of 



the jets and their appearance on the sky. The ephemeris used is listed in Table [TJ 



Our convention for labeling a piece of jet that we see at a certain spot on the sky is as 
follows: We define its birth epoch t as the time elapsed in the frame of the core since the 
material was born (ejected). Thus t serves to label specific pieces of the jets. We also define 
its age at emission r as the age of the same piece of material at the moment it emitted 
the photons that we detect at the same time as photons from the rest of the source. Due 
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to the finite speed of light, material that is "behind" the core and moving away from us 
had to emit its photons at a younger age in order that they arrive at the observer at the 
same time as photons from the core, and the reverse for matter "in front." This means that 
oncoming material has t > t, receding material the opposite. In this convention, the core 
has t = T = 0, material with larger t was born earlier than material with smaller t, material 
with larger r emitted its photons at an older age than did material with smaller r, and 
material moving in the plane of the sky has t = t. Quantitatively, birth epoch t and age at 
emission r are related by 

t 

l-v^/c 

where is the component of material's velocity that is directed toward the observer. In 
SS433, each jet has pieces with both positive and negative Vx, as illustrated in Figures [1] 
and El Figure [3] illustrates the light travel time effects in SS433. It shows the age at time 
of emission of photons r for each part of the east and west jets, as a function of the birth 
epoch t of each part of the jets. In the parts of the jets that are sufficiently bright that we 
can see them, the time differences {t — r) range up to ±200 days, grow with distance from 
the core, produce the well-known distortion of the appearance of the jets, and affect our 
analysis. Note that at most places along the jets, the differences between t and r for the 
two jets are not equal for a given value of t. 

8. APPENDIX 2 

Here we address the question of why the jets behave such that their observed 
brightnesses vary as D^~^°' rather than D'^^"', as might be expected for isolated individual 
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I 



components such as those seen in VLBI imagingo There are four issues that should not 
be confused. (1) Are the jets a series of isolated components or a continuous fluid flow? 
(2) For these two cases, how does the observed total intensity depend on Doppler boosting 
in a jet whose locus is a helix but whose motion is radial? (3) Under what conditions can 
observations distinguish the the cases n = 2 and n = 3? (4) What do the data say about 
the SS 433 jets? 

A moving optically-thin synchrotron source is Doppler boosted and K-corrected by a 



facto r of D^^°' because lu/i^^ is a Lorentz invariant and I^, oc u (iRybicki fc Lightman 



19791 ) . In the case of a continuous jet we want to know the flux densi t y of a segment of 



the je t defined by the o 



fll997h 



DeYoung 



^serving beam. As shown by 



Lind fc BlandfordI (119851 ) 



Sikora et al. 



( I2OO2I ), and others, not all the material in that segment of jet at a given 



instant contributes to a particular image. This is purely a light travel time effect. The 
fraction of the jet segment that contributes is (1 — (5 cos 6) = I/7-D. The effect of this is 
that the boost factor is modified from D^^"' to Z}^+"/7; in other words, J^*^* = 'jI°^^/D^~^"'. 
This calculation assumes a straight jet with the velocity along the jet direction. 

In SS433 the locus of the jet is a helix even though the motion of the jet material is 
radial. This means that the time delay between the arrival of photons from the near and 
far ends of the jet is no longer At = Rcos{6)/c, where R is the length of the part of the 
jet defined by the beam and 6 is the angle between the jet velocity and the line of sight. 
Instead, it is At ~ Rcos{ri)/c, where rj is the angle between the tangent to the jet locus and 
the line of sight. In addition, the rate at which material is added to the observed part of the 
jet is no longer proportional to Vjet, but is instead determined by the component of the jet 
velocity in the direction of the locus of the helix. We have made such calculations using the 



^The appearance of isolated components in VLBI images could simply be due limited 
dynamic range. 



- 15 - 



known kinematics of each part of the jet, and the results are very similar to the continuous 
straight jet case, that is, with effective boosts of approximately as observed. Critical 

to this analysis is the fact that the Lorentz factors of each part of the jets are the same. 

If the jets are each a series of components, then whether the total intensity is boosted 
by or depends on the spacing of the components relative to the spatial scale 

of the resolution of the imaging. If the components are sufficiently close, then there are 
always components that do not appear in the beam but belong there, and components that 
appear there but are not, just as for a continuous fluid flow, and Z}^+" is appropriate. If the 
spacing is so great that during the light travel time from the back of the jet to the front no 
component moves into or out of the beam, then is appropriate. However, depending 

on the properties of the components and the angular resolution of the observation, any 
value between n = 2 and n = 3 can be appropriate. For the current observation, the data 
strongly suggest that n = 2 for SS 433, indicating that the jets are either continuous, or if 
composed of discrete components, then there are many in the beam at any time. 
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Table 1. Ephemeris for SS433. 



Property 


SymboP 


Value 




R — 1) Ir 




1 1 t-L-CoolvJli 1 t-ilUU. 


P 




1— ; O T /_! T ' O T 1 1 ' t ^ H T ~\ / \ / ' M 1 1 1 1 1 


i 




Pypppcciori (/Orip Onpn i n o" A n o"1p 


lb 


20.92° 


Pypppccinn A^ic; TnpliriPifinn 


7 
L 


78 (15° 


rrecession Axis rosition Angle 




98.2 


Sense of Precession 


^rot 


-1 


Distance 


d 


5.5 kpc 


Orbital Period Reference Epoch (JD) 




2450023.62 


Orbital Period 




13.08211 d 


Nutation Reference Epoch (JD) 




2450000.94 


Nutation Period 




6.2877 d 


Nutation Amplitude 




0.009 


Porb 




0.0066 


Porb phase 




4.7 



'^As defined in 



Hjellming &: Johnston! (jl981bl ). 
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Fig. 1. — Total intensity images of SS433 made from 4.86 GHz VLA A-array data collected 
2003 July 11, using uniform weighting (ROBUST = —5). In (a) we show a contour image 
of SS433, with contours of total intensity /; the contours are spaced by factors of \/2; 
the peak is 301 mJy beam~^, the minimum contour is 50 /iJy beam~^, and the RMS 
noise is 21 /iJy beam~^. The kinematic model without jet velocity variations or nutation is 
shown as pluses for oncoming material and circles for retreating material, with components 
emitted at 5 day intervals, (b) Same as (a), except that the kinematic model including jet 
velocity variations and nutation is shown as blue and red lines, where these are oncoming 
and retreating parts, respectively. In these images the beam is 0.33" by 0.32". The scale of 
the figures is such that 1 " = 5500 AU, or 125 days of mean proper motion of 8 mas day~^. 
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Fig. 2. — Total intensity image of SS 433 made from 4.86 GHz VLA A-array data collected 
2003 July 11, using natural weighting (ROBUST = +5). We show contours of total intensity 
J; the contours are spaced by factors of a/2, the peak is 302 mJy beam~^, and the mini- 
mum contour is 26 /iJy beam~^, twice the root-mean-square noise of 13 /iJy beam~^. The 
kinematic model (not including jet velocity variations and nutation) is shown as blue and 
red lines, where these are oncoming and retreating parts, respectively. The beam is shown 
as a cross, and has full width at half maximum of 0.47". 
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Fig. 3. — Light travel time effects in SS433. The curves show the age of material in the jets 
as a function of the birth epoch of the material. The black line is for the east jet, the red 
line (grey in black & white) for the west jet. 
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Fig. 4. — Combined normalizations for the jets in SS 433, expressed as multiplicative factors 
to be applied to observed total intensities, as functions of the age of the jet material. These 
are the reciprocals of the beam-summed Doppler boosts for n = 2. The black lines are for 
the east jet, the reds lines (grey in black & white) for the west jet. Curves for n = 3 are 
similar and not shown. 
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Fig. 5. — Total intensity profiles of the jets of SS433 as a function of the age of the jet 
material derived from (a) the image in Figured! the root-mean-square noise is 21 /iJy/beam, 
and (b) derived from the image in Figure [21 the root-mean-square noise is 13 /xJy/beam. The 
black lines are for the east jet, the red ones (grey in black & white) for the west jet, and the 
broken lines show the core contributions to the total intensity. Jet emission is detectable out 
to ages at emission of at least 800 days (~ 6" from the core) in both jets. 
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Fig. 6. — Total intensity profiles normalized for (a) projection effects alone, and (b) for 
Doppler boosts {n = 2) alone (very similar to n = 3, not shown). The black lines are for the 
east jet, the red lines (grey in black & white) for the west jet, and the broken lines show the 
core contributions. 
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Fig. 7. — Intrinsic brightness profiles of tlie jets of SS 433 for {a) n = 2 and (b) for n = 3 as 
functions of the age of the jet material, derived from the image in Figure [H The black lines 
are for the east jet, the red lines (grey in black & white) are for the west jet, and the broken 
lines show the core contributions to the intrinsic brightness. 
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Fig. 8. — Observed total intensity ratio and normalization factor ratios, both in the form 
east jet divided by west jet, as functions of the age of the jet material r. The black line is the 
observed total intensity ratio, and the red line (grey in B&W) is the ratio of multiplicative 
normalization factors for n = 2 (Fig. Hj), the blue curve for n = 3, which provides a less 
satisfactory fit. If the two jets are identical and our normalization is correct, the red or blue 
curve should lie on top of the black curve. 
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Fig. 9. — Reconstruction of the west jet. We show the observed total intensities of east and 
west jets (black lines and red lines, respectively), and reconstructed total intensity profiles 
of the west jet (blue lines) as functions of the age of the jet material, for (a) n = 2 and (b) 
for n = 3. The data were obtained from the image in Figure [1] 
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Fig. 10. — Intrinsic brightness profiles of tlie jets of SS433 for (a) n = 2 and (b) for n = 3 
as functions of birth epoch, derived from the image in Figure [TJ The black lines are for the 
east jet, the red lines (grey in black & white) are for the west jet, and the broken lines show 
the core contributions to the brightness. 
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Fig. 11. — Linear, exponential, and power law fits to the intrinsic brightness profiles of the 
jets of SS433 as functions of the age of the jet material (for n = 2), for ages in the range 
60 < r < 150 days (data from Figured]), (a & b) Linear fits; the east and west jets have 
linear half-lives of about 46 and 40 days, respectively, (c & d) Exponential fits; the east and 
west jets have exponential half-lives of about 41 and 39 days, respectively, (e & f) Power law 
fits; the east and west jets have power law indies of about —1.7 and —1.8, respectively. In all 
parts the fits are the broken lines. As an example of the quality of the fits, the root-mean- 
square deviations from an exponential in the east and west jets are 0.8 and 1.0 mJy/beam, 
respectively. 



